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Phase diagram of the Hubbard-Holstein model in the coexistence of electron-electron and electron- 
phonon interactions has been theoretically obtained with the density-matrix renormalization group 
method for one- dimensional (ID) systems, where an improved warm-up (the recursive sweep) pro- 
cedure has enabled us to calculate various correlation functions. We have examined the cases of 
(i) the systems half-filled by electrons for the full parameter space spanned by the electron-electron 
and electron-phonon coupling constants and the phonon frequency, (ii) non-half-filled system, and 
(iii) trestle lattice. For (i), we have detected a region where both the charge and on-site pairing 
correlations decay with power-laws in real space, which suggests a metallic behavior. While pairing 
correlations are not dominant in (i), we have found that they become dominant as the system is 
doped in (ii), or as the electronic band structure is modified (with a broken electron- hole symmetry) 
in (iii) in certain parameter regions. 
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I. INTRODUCTION 

In condensed-matter physics, the electron-phonon in- 
teraction and electron-electron interaction are two of the 
fundamental interactions. While the problem is hard 
enough even when only one of them exists, it is a most 
challenging problem to consider what happens when both 
of the interactions coexist. 

A central problem is the ground-state phase diagram 
— how various phases arise in the coexistence of two 
kinds of interactions that arc quite different in na- 
ture. Of particular interest is superconductivity. Strong 
electron-electron repulsive interactions introduce spin 
fluctuations in the electronic system, which can medi- 
ate pairing interactions for the electrons to make them 
condense into superconducting states with anisotropic 
Cooper pairs?i The superconducting phase has to, how- 
ever, compete with density- wave phases such as the spin- 
density wave (SDW), which also come from the electron- 
electron interaction.^ On the other hand, the electron- 
phonon interaction mediates an attraction between elec- 
trons, which can make the electrons condense into super- 
conducting states with isotropic Cooper pairs. This time 
the superconducting phase has to compete with density- 
wave phases such as the charge-density wave (CDW) aris- 
ing from the electron-phonon interaction. The electron- 
phonon interaction is the coupling of the conduction elec- 
trons to the lattice structure, and it can also work to 
make the system undergo a Peierls transition, where the 
lattice is deformed so that the electrons are not conduct- 
ing. 

So it is quite a nontrivial problem to consider the 
ground-state phases when the two interactions coexist 
and are both strong. A prototypical model represent- 
ing such a situation is the Hubbard-Holstein model, 
where the Hubbard model for the electron correla- 
tion is coupled to (Einstein) phonons. The model is 



characterized by three physical parameters: (a) the 
on-site electron-electron repulsion U , (b) the phonon 
frequency wq, and (c) the electron-phonon coupling 
A, where the unit of energy is the electronic transfer 
t cx electronic bandwidth. Of particular interest is 
the intermediate regimes: (a) the regime {Tiloq ^ t) 
intermediate between the adiabatic {Tiuoq <C W) and 
anti-adiabatic [Huiq :$> W) limits, (b) the regime {U ~ A) 
where the electron-electron Coulomb repulsion and 
the phonon-mediated attraction are similar in magni- 
tude. The problem has in fact been studied in various 

approaches ,3-4.5. 6. 7.8. 9 . 10.11.12. 13.14. is. 16.17. is. 19.20. 21. 22. 23. 24. 25. 26 

This is by no means a theoretical curiosity, since 
there are various classes of materials in which both the 
electron-electron and electron-phonon interactions are si- 
multaneously strong. A typical example is the solid 
fuUerene doped with alkali-metal atoms, (A3 Ceo, ^ — 
K, Rb, where the Ceo fuUerene molecules are 

aligned in a face-centered cubic (fee) lattice, for which 
the alkali atoms supply electrons making the conduc- 
tion band half-filled. In this material the intramolec- 
ular phonon modes are known to have high frequencies 
0.2 eV), which couple strongly to the conduction elec- 
trons. The material is a superconductor whose transition 
temperature Tc is the highest to date among the car- 
bon compounds. It is an interesting problem to consider 
how the half-filled electron band becomes metallic, rather 
than becoming CDW or SDW insulators from coexisting 
electron-electron and electron-phonon interactions, and 
exhibits superconductivity. The electronic bandwidth of 
A^Cqo can be controlled by e.g. changing the alkali metal 
species A, where Tc is observed to increase with the 
lattice constant. Organic chemists have also fabricated 
fuUerene derivatives in the shape of a shuttlecock, where 
the molecules stack to form an array of ID chains?^ 
where the separation between Ceo molecules in a chain 
is controllable with the choice of functional groups at- 
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tached to Cgo, and the system is shown to have stronger 
electron correlation,^- 

Given these backgrounds, the purpose of the present 
paper is to obtain the ground-state phase diagram when 
the electron-phonon and electron-electron interactions 
coexist over the whole parameter space spanned by 
U, X, huj that includes the intermediate regimes (Hlu ^ W] 
U ^ X). We adopt the Hubbard-Holstein model, one of 
the simplest models for the coexisting electron-electron 
and electron-phonon interactions. We want to treat these 
interactions on an equal footing, especially in the inter- 
mediate regimes. Phonons mediate attraction between 
electrons, but here we have adopted the dispersionless, 
Einstein phonons that do not directly propagate along 
the chain. While these phonons cannot directly mediate 
the pairing force between electrons across different sites, 
we can still expect the phonons to contribute to the pair- 
ing, because electrons hop between neighboring sites. So 
the phonon self-energy should depend not only 

on Lo but also on q. Specifically, phonons with q ^ n/a 
(a: the lattice constant) as well as phonons with q ^ 
should be strongly renormalized by the coupling to the 
electrons, because the Fermi points reside at fc = ±7r/2a. 
This is another reason why the electron and phonon de- 
grees of freedom should be treated on an equal footing. 

So a special care has to be taken in choosing the 
method, since the method has to (a) take account of 
the phonons without assuming the adiabatic or anti- 
adiabatic limits, and (b) take account of the electrons 
where charge gaps around the half-filled electronic band 
(a region of interest) can be described. We have adopted 
the density-matrix renormalization group (DMRG) as an 
appropriate method. 

DMRG was originally developed for spin systems22ii2^ 
and applications to electron systems subsequently 
followed. ■^"''■^^ The basic idea is the following: We iter- 
atively add sites to enlarge a one-dimensional system, 
which is approximately represented as two connected sub- 
chains with truncated Hilbert spaces optimized for a 
target state of the whole system. In adding a site to 
a subchain, we calculate the partial density matrix for 
the target state wave function, diagonalize it, and re- 
tain eigenstates that correspond to a limited number of 
largest eigenvalues, to represent the new Hilbert space 
for the new subchain. While the DMRG was developed 
for electron systems, Jeckelmann and White extended 
the formalism to incorporate phonons in DMRGj^^ We 
have here adopted this formalism, with an improved al- 
gorithm. 

In constructing a phase diagram, a most direct method 
is to actually calculate various correlation functions. 
However, the inclusion of phonons, which are bosons, 
makes the calculation enormously heavy, and we have to 
somehow overcome this difficulty. Here we have adopted 
an improved warm-up (the recursive sweep) procedure for 
the DMRG developed by one of the present authorSf^ 
which has enabled us to actually calculate correlation 
functions (including the pairing correlation) for the first 



time for the Hubbard-Holstein model. While taking ID 
systems is a theoretical device (apart from the above 
mentioned quasi- ID fuUerides), it is theoretically intrigu- 
ing to look at functional forms of correlation functions in 
the coexistence of electron-electron and electron-phonon 
interactions. The Mermin- Wagner theorem dictates that 
continuous symmetries of the Haniiltonian are not broken 
in the ground states of infinite ID systems. Hence the 
correlation functions for the phases with broken contin- 
uous symmetry should decay with distance. We can still 
compare the exponents of the decay, as in the Tomonaga- 
Luttinger theory for purely electronic systems, to deter- 
mine the dominant correlation, which may become long- 
ranged once we go to higher dimensions. 

Now, there are many existing theoretical studies on the 
Hubbard-Holstein model, but the nature of the model has 
yet to be fully established. There have been competing 
results even on the ground-state phase diagram at half- 
filling. Particular issues are: 

1. For large enough U we can expect the ground state 
to be an SDW, while for large enough A we can ex- 
pect the ground state to be a CDW. An interesting 
problem is whether there exists a metallic phase 
between the two density-wave phases. 

2. Whether a nai've expectation is valid, where one 
imagines that, for the coexisting electron-electron 
and electron-phonon interactions, the net interac- 
tion will simply be the difference between them. 
This problem becomes especially interesting in the 
intermediate regime. 

3. Whether the metallic phase, if any, can possibily 
become superconducting. For pairing symmetries, 
we have here considered all the possibilities of 

• sSC: on-site spin-singlet pairing, 

• pSC: nearest-neighbor spin-triplet pairing, 
and 

• dSC: nearest-neighbor spin-singlet pairing. 

It may sound peculiar when we say p- or d-wave 
superconductivity for ID systems, but a nearest- 
neighbor spin-triplet (singlet) pairing corresponds 
to a p((i)-wave pairing in two or higher dimensions, 
so we adopt this nomenclature. 

4. How the phase diagram depends on the electron 
band filling, or the lattice structure (i.e., the elec- 
tronic band structure) of quasi-lD systems. 

In order to resolve these issues we have obtained the 
correlation functions for long chains with DMRG for the 
cases of (i) the systems half-filled by electrons for the full 
parameter space spanned by the electron-electron and 
electron-phonon coupling constants and the phonon fre- 
quency, (ii) non-half-filled system, and (iii) trestle lattice. 
For (i), we have detected a region (where the electron- 
phonon interaction is stronger than U but not too strong) 
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where both the charge and on-site (sSC) pairing correla- 
tions decay arithmetically in real space, which suggests 
a metallic behavior. While pairing correlations are not 
dominant in (i), we have found that they become domi- 
nant as the system is doped in (ii), or as the electronic 
band structure is modified (with a broken electron-hole 
symmetry) in (iii) in certain parameter regions. 

Organization of the paper is as follows: In section II 
we introduce the Hubbard-Holstein model, and briefly re- 
view the open questions about the phase diagram. In sec- 
tion III we explain the numerical method for calculating 
correlation functions. Results are presented in sections 
IV (half-filled chain), V (doped chain) and VI (undoped 
trestle lattice). Concluding remarks are given in section 
VII. 



II. MODEL 




-t 



FIG. 1: Tlie Hubbard-Holstein model. Each atom consists of 
an electron site (circle) and a local phonon with frequency ujq 
(spring). Arrows denote up and down spin electrons, which 
are coupled to the phonons with a strength A. An electron can 
hop between neighboring sites with the hopping amplitude t, 
and feels an on-site electron-electron repulsion when encoun- 
tered with another electron of opposite spin on the same site. 



We take the Hubbard-Holstein model, where the elec- 
tronic part of the Hamiltonian is the Hubbard model that 
in itself exhibits various phases according to the values 
of the electron band filling n and the short-range (on- 
site) electron-electron interaction U, dimensionality, and 
the lattice structure. Phonons are introduced as a local 
harmonic oscillator for each site, to which the electrons 
are coupled as shown schematically in Fig. [1] Thus the 
phonons are dispersionless (i.e., Einstein phonons). In- 
clusion of on-site phonons adds two parameters to the 
Hubbard model: the electron-phonon coupling g and the 
phonon frequency loq. We can characterize the strength 
of the electron-phonon coupling, 

A = 2g'^/huJo, 

which is the attraction mediated by phonons between 
two electrons on the same site in the anti-adiabatic 
{u!o/t oo) limit. We take this as the measure of the 
electron-phonon coupling to facitilate comparison with 
the electron-electron repulsion U. We set fi = 1 here- 
after. 

The Hamiltonian is given as 

+ gniajai + a|) + ^a;oa|'a,. (1) 

z.cr i 

Here, Cia- annihilates an electron with spin (t(=|,4) at 
site i, ni„ — c\^Ciu is the electron number, and hi is the 
phonon annihilator at site i. If we consider the phonon 
as a harmonic oscillator where a mass M is attached to 
a spring with a spring constant K with the phonon fre- 
quency uoq = K/M ^ the displacement Xi is written as 



The electron-phonon coupling is introduced as a coupling 
of the local electron number, Ui = X^o-"*'^' with the lat- 
tice displacement Xi. 

The electron band filling is definded a&n„ ~ and, 
by introducing a = a — gn„ /Tiloq, we can rewrite ([T]), up 
to a constant, as 

i 

+ ^^g{ni„ -na){ai + a\) + ^^hujaa\ai. (3) 

z.cr i 

In this form, the relevant electron occupation is the de- 
viation from the average value. 

The Hubbard-Holstein model ([T]) has three indepen- 
dent parameters: 

U/t,\/t,uja/t 

as schematically displayed in Fig. [21 In this parameter 
space, we are interested in the regimes (i) away from the 
adiabatic (o^o/i ^ 1) or anti-adiabatic {oj^/t ^ 1) limits, 
and (ii) comparable U ^ X. To summarize the existing 
theoretical results on the ID half-filled Hubbard-Holstein 
model, we observe the following issues; 

• Two of the recent studies on this model, one with 
the Lang-Firsov transformation'^^ and the other 
with DMRG,^^ have different conclusions as to 
where this model becomes metallic. The former in- 
dicates a considerably wide metallic region between 
the CDW {U < A) and SDW (A < U) phases, 
while the latter indicates a closing of the charge 
gap only at a quantum critical point between the 
CDW and SDW insulators. 

• While the strong-coupling expansion^ indicates 
that superconductivity cannot be dominant for fi- 
nite ujQ, a study with some ansatz for phonons^ 
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indicates a dominant superconductivity near the 
U — axis on the U-X phase diagram, and an- 
other QMC study for the charge structure factor 
suggests a phase diagram with a superconducting 
region near the U = X hne between the CDW and 
SDW phases. 

On the other hand, the Holstein model corresponds to 
the [/ = plane in the phase diagram of the Hubbard- 
Holstein model, and picks up the effect of the electron- 
phonon interaction in making the electronic system insu- 
lating. In one^i^ and two dimensions, and in the limit 
of infinite dimensionj^i^ there are extensive studies on 
the transition between (i) the small-polaron regime where 
the electrons are self-trapped by the lattice through a 
stronger electron-phonon interaction, and (ii) the large- 
polaron regime where the electrons move around ex- 
citing local phonons through a weak electron-phonon 
interaction!^ When the electron-electron interaction is 
turned on in such an electron-phonon system, the in- 
terplay between the effects of the electron-phonon and 
electron-electron interactions again becomes interesting, 
which has not been fully understood. 

III. METHOD 

As we have seen, the ground-state phase diagram for 
the ID Hubbard-Holstein model is controversial. The 
most clear-cut way for examining the competition of var- 
ious phases is to look at the correlation functions on long 
systems, which has not previously been done properly. 
So here we calculate correlation functions in real space. 
A care has to be taken in examining correlation func- 
tions in ID quantum systems: strong quantum fluctua- 
tions repress continuous symmetry-broken long-range or- 
ders in the ground state.— However, if one type of 
correlation function decays more slowly than others, we 
can identify the correlation to be dominant. Exact re- 
sults for purely electronic systems on various integrable 
models show that diagonal [{c^ {r)c{r)c^ {Q)c{Q))] and pair 
[(c^(r)ct(r)c(0)c(0))] correlations decay with a power-law 
(r"'', where rj is an exponent) or exponentially (e"''/^, 
where ^ is the correlation length) at large distances r. 

In the case of charge and spin correlations, the charge 
(spin) correlation function decays with a power-law when 
the charge (spin) excitation is gapless, Achargc — 
(Aspin = 0), while the correlation decays exponentially 
when the charge (spin) excitation has a gap Achargc > 
(Aspin > 0). To be precise, these correlations can have 
logarithmic corrections or additional terms that decay 
faster than the main term. This is one reason why we 
need long chains to identify the dominant phase by nu- 
merically calculating the decay of correlation functions. 
Here we compare the real-space behavior of various cor- 
relation functions in the ground state for a chain having 
as many (typically 64) sites as numerical calculation is 
possible for model parameters away from adiabatic and 
anti-adiabatic limits. For these reasons, we adopt the 



DMRG with the pseudo-site method^ to calculate the 
ground-state correlations. We adopt the open-boundary 
condition for a better DMRG convergence. The recursive 
sweep initialization'^^ adds two sites for both the electron 
and phonon degrees of freedom at each infinite-algorithm 
DMRG step, while benefiting from the reduced size of the 
Hilbert space in the pseudo-site method. This also con- 
tributes to better convergence with less computational 
resources, which has allowed us to compute correlation 
functions on lattices with large number of sites. 



IV. RESULT: HALF-FILLED 
HUBBARD-HOLSTEIN CHAIN 

A. Parameter space 

First we investigate half-filled systems for the param- 
eter space shown as a bird's eye view in Fig. [2] (with 
t — 1 taken as the unit of energy), while the doped case 
is treated in f|Vl 

We calculate charge, spin and pair correlation func- 
tions for this parameter space. For the pairing symmetry 
we consider the on-site, spin-singlet pair (sSC) (which 
corresponds to s-wave superconductivity in higher spa- 
tial dimensions) , spin-singlet pair across the neighboring 
sites (dSC) (d-wave in higher dimensions), and triplet 
pair across neighboring sites (pSC) (p-wave in higher di- 
mensions). Their operator forms are: 

• Charge: {'n^nj) - {nt){nj), 

• Spin: (SfSj), 

• sSC: i^l^j), 
. pSC: (hIh,), 

• dSC: (A,|Aj), 

where = = (n^i - ni{)/2, = Ci^Cix, 

Hi = (ci-fCj+ix -I- Ci^Ci+i^)/\/2, and Ai = {a^Ci+ii - 

C^iC^+^)/V2. 

We first display typical dependence of correlation func- 
tions on the real-space distance along the chain. The 
calculations have been done on L = 64-site, half- filled 
Hubbard-Holstein chains, with pseudo-site DMRG for at 
least Nf, = 3 phonon pseudo-sites per full site, unless oth- 
erwise indicated. At least 10 finite-system sweeps have 
been done with m > 400 states retained in more than 
two final iterations after the infinite-finite method warm- 
up. The results do not change significantly when we take 
larger Nb > 4 or m > 500. The maximum discarded 
weight of the partial density matrix in the final sweep 
is below 10~^ when m > 400 and typically around 10~^ 
when m > 500. 
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FIG. 2: Parameter space for the half-filled Hubbard-Holstein 
model spanned by the on-site electron-electron repulsion U , 
the phonon frequency tJo, and the electron-phonon coupling 
A, where the unit of energy is the electronic transfer t. The 
model reduces to the Hubbard model when A = 0, and to 
the Holstein model when U = Q. The uio/t limit is the 
adiabatic limit, while the o^o/i — > oo limit is the anti-adiabatic 
limit. 



{X, coq) = (3.6, 5.0) 
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FIG. 3: (Color online) Log-log plot of correlation functions 
against the real-space distance for (A,tJo) ~ (3.6,5) and U = 
(a) 8.0, (b) 3.6 and (c) 1.6 for the half-filled Hubbard-Holstein 
model. 



B. The case of larger cjo/t 

First we consider the region where the phonon fre- 
quency is greater than the electron hopping t. We fix 
the electron-phonon couphng A and phonon frequency 
ujo at (A,(jJo) = (3.6,5.0). We observe in Fig. [Sja) that 
only the spin correlation decays with a power-law for U 
much larger than A. The exponent is near unity. 

As U is decreased to around A, the charge correlation 
begins to decay with a power-law, as plotted in Fig. E^b) 
for J7 = 3.6 = A. The exponent for the spin correlation is 
greater than in Fig. [21(a). Here we note that the on-site 
pair correlation also decays with a power-law. 

Finally, for a much smaller value oi U = 2.6 in Fig. 
^c), the spin correlation decays no longer with a power- 
law, but almost exponentially at large distances, which 
becomes more manifest as U is further decreased. There, 
both CDW and on-site pair correlations still decay with 
power-laws. So we observe a power-law spin correlation 
when the electron-electron interaction is stronger than 
the electron-phonon interaction, or power-law charge and 
on-site pair correlations when the electron-electron inter- 
action is weaker. 

Next we compare the exponents of the correlations. We 
fit the correlations as functions of the real-space distance 
r, assuming the behavior {0^{x)0{x + r)) ^ r"'', and 
determine the exponent rj. The results are plotted in Fig. 
131 We can see that charge and on-site pair correlations 
have similar exponents, but the exponent for the charge 
correlation is always smaller than that for the on-site pair 
correlation. So the charge correlation is dominant when 
U/t is smaller than around 3.5 where A = 3.6, while the 
spin correlation is dominant when U/t is larger. 

For the repulsive {U > 0) Hubbard model, we have 
only the spin correlation decaying as r~^. Now, for the 
purely electronic Hubbard model, it has been known^i^ 
that we can convert the repulsive {U > 0) Hubbard 
model into the attractive {U < 0) Hubbard model, at 
half-filling, with a canonical transformation, with which 
the charge and sSC pair correlations are mapped to the 
spin correlation in the repulsive side as 

U ^ -U 

« {srs+). (4) 

This implies that the charge and sSC pair correlations are 
degenerate and both decay like as far as the Hubbard 
model {U < 0) and the Hubbard-Holstein model {U = 
[/ — A < 0) in the anti-adiabatic limit are concerned. 

If we go back to the Hubbard-Holstein model, the 
power-law fit for an exponentially decaying correlation 
gives a finite exponent for finite systems, so what we 
should observe when a correlation does become power- 
law is a decrease of the exponent around the phase tran- 
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{X, coq) = (3.6, 5.0) 
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FIG. 4: (Color online) Exponents of correlation functions for 
the Hubbard-Holstein model with (A,tJo) ~ (3.6, 5.0). 
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FIG. 5: (Color online) Phonon number {nphonon(i)) plotted 
against i for U — 1.8,3.6,4.2 and {X,ujo) — (3.6,5.0) for a 
64-site, half-filled Hubbard-Holstein chain. 



Phonon numbers per site in the 
Hubbard-Holstein model 



We can calculate the phonon occupation number for 
each site nphonon(«) in the ground state Tpo, defined as 



(5) 



We can check the boundary effects by plotting the 
phonon number against the site number i in Fig. [5l The 
curves are almost flat except for a few sites around either 
edge of the chain. 

We plot the phonon number at the center of the chain, 
('T-phonon), sls & function of U in Fig. ^i). Phonons me- 
diate an attraction between two electrons with opposite 
spins on the same site, which amounts to A = 2g^ /ujq in 
the loq ~f oo limit. For a finite loq, the attraction A is 
smaller than this asymptote A,"*^'^^ For U X, two elec- 
trons tend to form a bipolaron, where the electrons are 
strongly bound to each other by sharing phonons they 
excite. Bipolarons, at half-filling, tend to occupy every 
second site to reduce kinetic energy. In this case the 
charge correlation that corresponds to the 2kF CDW is 
the strongest. The number of excited phonons per site 
increases as U is decreased, because the fluctuation from 
the single occupancy increases as the CDW correlation 
becomes stronger. 

The limit of immobile bipolarons can be understood 
as follows: we consider two electrons on a site in ([3]) 
and neglect their hopping, to apply a Lang-Firsov type 
transformation d ~ a — y/a. Then the local Hamiltonian 
for the phonon on this site is 

H = g{a + a^) + LUoo'a = ujoo^a + {ujQa + 2g\/a), (6) 

with a = (g/ujQ)'^ = X/2ujq. The second term is a 
constant, and if we denote the n-phonon state as \n) 
(n = 0, 1, . . .), the ground state IV^o) = ^"1'^) satisfies 

a\ipo) = {a~y/a)'^bn\n) = y^(\/n6n-\/a)|n-l) = 0. 

n n>l 

(7) 

This is solved as a coherent state, 



sition point. We observe in Fig. 2] that the exponent 
for the SDW (CDW and sSC) decreases (increases) as 
[/ = [/ — A is increased from the negative to positive 
sides. This behavior is similar to the Hubbard model 
having U as the on-site interaction, except that the ex- 
ponent for the charge correlation is slightly smaller than 
that of the sSC pair correlation in general for finite cuo/t 
as mentioned above. This is similar to the case of the 
Holstein model U = with an LUo/t ^ 5 (not shown) . 



which gives 



lim (nphonon) = (a^a) = 



a = \/2ujQ. 



(8) 



(9) 



We can observe in Fig. WLi) that the system indeed ap- 
proaches to this limit as U is decreased. 

A deviation observed for larger values of wq in Fig. [6Uf ) 
is interesting in that the region is intermediate between 
the J7 <C A and U ^ \ limits. There we can expect a 
behavior different from the limiting behavior. 
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X = 3.6 



(U, X, coq) = (2.1, 3.6, 2.0) 
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FIG. 6; (Color online) (a) - (e) Exponents of correlation 
functions determined by calculating exponents for a 20-site 
Hubbard-Holstein model for \/t — 3.6 with various values of 
ujo/t = 0.05 — 5. (f) The phonon number per site, normalized 
by eqn.®. 
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FIG. 7; (Color online) Correlation functions for {U,\,u)) = 
(2.1,3.6,2.0) on a 64-site, half-filled Hubbard-Ifolstein chain 
with ub ~ 6, where CDW correlation is dominant and almost 
constant. 



D. The case of smaller uJo/t 

When the value of u!o/t is smaller than about 3, we 
start to observe in Fig. [6fc) behaviors of the correlation 
functions different from those discussed above: the value 
of U at which the exponent for the charge and spin corre- 
lations coincide becomes considerably smaller than A. So 
the (U—X) Hubbard model picture is no longer valid. Sec- 
ondly, the charge correlation becomes almost flat against 
distance as in Fig. [71 implying a charge-ordered (CO) 
phase, in some region of U in Fig. [nila)-(d), while this 
does not happen for the on-site pair correlation. 

To analyze the nature of this behavior, we can again 
look at the phonon occupation number against U for 
various values of ujq for a fixed value of A (Fig. HUf)). 
For small enough loq the phonon occupation number 
abruptly increases at some value Uc as U is decreased. 
We can see that the exponential decay of all the correla- 
tions except the charge correlation in fact start around 
U = Uc- Thus we identify that this is where the SDW 
- CDW transition occurs. For larger uiq, on the other 
hand, the change in both the exponents and the phonon 
occupation number are gradual. The intermediate val- 
ues (0 < (jT-phonon) < A/2ti;o) of the phonon occupation 
number is consistent with the power-law behavior of the 
CDW and sSC exponents, which corresponds to a metal- 
lic phase with a closed charge gap, where the fluctuation 
in the local electron number (hence the fluctuation in the 
local phonon number) will be large. 

When we look at the CDW correlation. 



Ccharge(i, j) 




^0 



(10) 

its amplitude turns out to either decay, or remain al- 
most constant (^co) against the distance r = \i — j\. 
Since {J2a '^i") alternates between 1 -I- /co and 1 — /co 
in the latter case at half-filling, we can simply define 
the charge-ordering amplitude /co = \/|^co|- We plot 
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{X, coq) = (3.6, 2.0) 



CD 0.8 

T3 

i. 0.6 

g 0.4 
D 
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FIG. 8: (Color online) Charge gap, spin gap, and CDW am- 
plitude plotted against U for the half-filled Hubbard-Holstein 
model with (A, wo) = (3.6,2.0). Data obtained for L = 
20, 30, 40, 64 have been linearly extrapolated to 1/L ^ 0. 



Achargc, Aspin and fco against U in Fig. [S] As C/ is de- 
creased in the CO region, we observe that (i) the CDW 
amplitude /co increases and approaches unity, (ii) both 
the spin and charge gaps increase almost Unearly, and 
(iii) the spin gap is nearly twice the charge gap. For the 
CDW/sSC region, the charge gap almost vanishes. The 
spin gap becomes also small, but remains larger than the 
charge gap. Finally, when U exceeds A, the spin gap 
remains zero, while the charge gap slightly grows with 
U. 

The effect of the ratio ujo/t on the ^/-dependence of 
(riphonon) can now be understood in the limits of S> t 
and ojQ ^ i as follows: 

1. In the anti-adiabatic limit (wo ^ t), the response 
of the phonons to the motion of electrons is so fast 
that we can think of the role of phonons as just me- 
diating the attraction. A, between electrons. Elec- 
trons then tend to occupy different sites for [/ > A, 
with spins antiparallel between neighbors because 
of the exchange interaction. In this case the spin 
correlation that corresponds to the 2kp SDW is the 
strongest, and the number of excited phonons per 
site decreases with U, because the electron occu- 
pancy becomes closer to unity for every site. 

2. In the adiabatic limit {ujq ^ i), for t/ <C A the de- 
viation (/co) of the electron occupancy from the 
average value is large, so that a larger number of 
phonons are excited. The dynamics of phonon is 
much slower than that of electrons for uq t, 
and exactly when the charge-ordered states exist 
depends on as we shall show on the phase dia- 
gram. For U ^ X electrons cannot move and form 
a Mott insulator, so that phonons are hardly ex- 
cited. In the region loq ^< t we observe a jump 
in (riphonon) as a function of U . This should corre- 
spond to the CDW transition. 



TABLE I: Summary of possible phases at half-filling. We 
denote power-law and exponential decays of correlation func- 
tions ('corr.') as 'power' and 'exp.', respectively. 



phase 



CO 



CDW/sSC SDW 



spm corr. exp. exp. power 

charge corr. long-range order power (smallest 77) exp. 

sSC corr. exp. power exp. 

charge gap finite finite 

spin gap finite finite 

(wphonon) — > \/2uJo intermediate ^ 



E. The phase diagram 

As we have seen above, we have two to three distin- 
guishable regions out of the five phases compared in Table 
Uwhen we increase the value of U from [/ = 0: 



1. Charge ordering (CO) (A 



spin 



A, 



charge 



> 0) 



The charge correlation is dominant and does not 
exhibit a significant decay against distance. Other 
correlation functions decay exponentially. The 
charge and spin gaps are still large, but decrease 
with U . These should be the characteristics of 
a long-range charge order, which is not forbidden 



frc 



^43, 44 



We denote this region as CO. 



2. CDW/sSC (Achargc < i,0 < A,pi„ < t) — The 
charge and sSC correlations decay with power-law, 
with the charge correlation the stronger of the two. 
The spin gap is finite (although a size-scaling argu- 
ment will be required to quantify this), while the 
charge gap is closed (although a size-scaling argu- 
ment will again be required). We denote this region 
as CDW/sSC because the sSC correlation, while 
not dominant anywhere, decays with power-law as 
well. The dSC correlation (nearest-neighbor spin- 
singlet pair) is even less dominant. 

3. SDW (Achargc > 0, Aspin <C i) — Only the spin cor- 
relation decays with power-law with the exponent 
being nearly unity. The spin gap is closed or much 
smaller than t. The charge gap becomes larger as 
U is increased. We denote this region as SDW. 

We draw phase diagrams on the t/ — A plane for various 
values of = 0.2,0.5,2, and on the U ~ uia plane for 
various values of A = 0.65, 1.5, 3.6, 8.4 (Fig. [Eb)(c)). We 
observe on the U — luq plane that the CDW/SC region 
is narrow for ^ A, while the region expands with the 
CO-CDW/SC boundary shifting to smaller U and the 
CDW/SC-SDW boundary shifting to larger [/, until the 
region finally extends up to f7 ~ A when loq >~ A. The 
CO region also becomes wider for smaller loq or smaller 
U, but vanishes at ojq ~ A. For A — 0.05 we have not 
detected this region. 
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CDW/SC 




(U, X, coq) = (2.0, 3.6, 5.0) 



° \CDW/SC U 



FIG. 9: (Color online) (a) The summary phase diagram on the 
whole parameter space for the half-filled Hubbard-Holstein 
model as derived from (b) the results obtained for the U — X 
cross sections, and (c) U — logojo cross sections. The bound- 
aries have been determined from the data points marked with 
dots. 



On the U — \ plane, we observe that the SDW region 
covers all oi U > A, but it extends beyond the line U = X 
when ojQ is smaller. The CDW/sSC region is broader for 
larger values of wq, and we have not detected CO region 
for LUo/t = 20. Figure (nja) summarizes the whole region. 
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FIG. 10; (Color online) Correlation functions for a 64-site, 
hole-doped [n = 30/64 = 0.47 in (a) and n = 28/64 = 0.44 in 
(b)] Hubbard-Holstein chain with {U,X,loo) = (2.0,3.6,5.0). 
The exponents for the sSC correlation (dotted straight line) 
are estimated to be t^sSC ~ 1-16 ± 0.02 in (a) and ?7ssc = 
1.09 ± 0.02 in (b), while the CDW correlation to be »7cdw = 
1.30 ± 0.10 in (a) and r?cDW = 1-54 ± 0.13 in (b). 



V. DOPED CHAIN 

Having seen that the pairing correlation is only sub- 
dominant in the half-filled Hubbard-Holstein model, we 
can naturally ask the question: can the system away from 
the half-filling be superconducting? Since the effect of 
doping is dramatic in purely electronic systems such as 
the Hubbard model, this is an obviously interesting av- 
enue to examine. For the doped systems we change in our 
DMRG calculation the target quantum number, which 
is the pair of the up spin and the down spin electrons 
(7V|, Ni), and change the values of rTo- in ^ accordingly. 

The functional forms of the correlations for the non- 
half-filled system have turn out to have the following 
property: the pair correlation functions decay almost 
exactly with a power-law, while the charge correlation 
decays as a power-law multiplied by an oscillating fac- 
tor. The dominant wavenumber for the CDW correla- 
tion should be 2fci?, which equals to ir/a (a: the lattice 
constant) at half-filling, but becomes larger as the hole 
doping is increased. In fact, we can fit the result for the 
CDW correlation as 



CcDW ^ Cr- ^c^* cos {qr + c) 



(11) 



where TycDW is the exponent, C, c and (7(~ mr/a) are 
fitting parameters. 

We compare the CDW exponent thus obtained with 
the exponent for the sSC correlation (dashed line) in Fig. 
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Trestle, (U, X, coq,?') = (2.5, 3.6, 5.0, -0.4), n = 1 




r 



FIG. 11: A trestle lattice (bottom) is equivalent to a chain 
that has the next-nearest neighbor hopping t' (top). 



[TUl The calculated exponents for CDW and sSC are sim- 
ilar in Fig. [10] (a), where the electron filling n is close 
to unity. However, sSC becomes clearly dominant over 
CDW for a larger level of hole doping {n = 0.44) in Fig. 
[Tn](b). The exponents in this case are: ?7ssc — 1.09±0.02 
and 77CDW = 1.54 ± 0.13. So this is the key result for the 
doped system: the sSC correlation has a smaller expo- 
nent, and thus dominant in the sufficiently doped case. 

The doped system further poses an interesting prob- 
lem: for purely electronic systems such as the Hubbard 
model, doping the Mott insulator in two or higher dimen- 
sions favors superconductivity. For the Hubbard-Holstein 
model with phonons, however, the situation should be far 
from trivial, especially because the simplified U = U ^ X 
picture can be invalidated as stressed above. So this 
should be a further avenue for future studies. 



VI. EFFECT OF THE LATTICE STRUCTURE 
— TRESTLE LATTICE 

A second approach to make the pairing correlation 
dominant should be, in our view, to modify the electronic 
band structure. Along the line of argument above, we can 
specifically break the electron-hole symmetry to lift the 
degeneracy between CDW and sSC correlations, eqn.([5]), 
which exists in the Hubbard modelM and is shown above 
to persist for the Hubbard-Holstein model in the anti- 
adiabatic limit. This can be done by considering quasi- 
ID lattices such as ladder or trestle lattices. A trestle 
lattice, for instance, is equivalent to a chain where a next- 
nearest neighbor hopping t' is introduced (Fig. [TT|) . This 
washes out the bipartite symmetry in the band structure, 
which in turn breaks the electron-hole symmetry. We can 
then seek the possibility of making the pairing correlation 
dominant, even at half-filling, in the Hubbard-Holstein 
model, as have been done for the Hubbard model4^'^ 

So we have performed the DMRG study for the tres- 
tle lattice. We retained up to m = 720 states per 
block. The maximum truncation error in the final sweep 
is around 10^"", so it is considerably larger than in the 
case of the simple chain. If we look at the correlation 
functions against the real-space distance in Fig. [121 we 
can see that the sSC correlation is indeed dominant at 
{U, X,ujQ,t') = (2.5,3.6,5.0,-0.4), even though we con- 
sider a half-filled system. Here both of the charge and 
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FIG. 12: (Color online) Correlation functions plotted against 
the real-space distance r for a 64-site, half-filled Hubbard- 
Holstein model on the trestle lattice with {U, X,uJo,t') = 
(2.5,3.6,5.0,-0.4). 



Trestle, (A., coq) = (3.6, 5.0), n=1 
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FIG. 13: (Color online) Exponents of correlation functions 
plotted against U for a 64-site, half-filled Hubbard-Holstein 
model on the trestle lattice for {\,uio) — (3.6, 5.0) with t' /t — 
—0.3 (top panel) or —0.6 (bottom). 



sSC correlation functions exhibit power-law decay times 
oscillating factors. We can then look at the exponents of 
correlation functions against U for various values of t' /t, 
which controls the band structure (top insets of FigfT4|). 
A special interest is that the number of Fermi points at 
half-filling for the non-interacting system increases from 
two to four at \t'\ > 0.5. Curiously enough, the case of 
four Fermi points is seen in the bottom panel of Fig. [T3] 
to have a dominant dSC correlation for U X. This is 
the key results for the trestle lattice. 

For U X we observe relatively larger error bars in 
the exponents. This should be due to the presence of 
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Trestle, {U, X, coq) = (2.0, 3.6, 5.0), n = 1 




FIG. 14: (Color online) Correlation exponents plotted against 
—t' for a 40-site, half-filled Hubbard-Holstein model on the 
trestle lattice with {U,X,iJo) — (2.0,3.6,5.0). At least m — 
630 states have been retained in the last (10th) sweep of the 
finite algorithm DMRG. The vertical dashed line represents 
the boundary at which the number of Fermi points changes 
from two to four as indicated by the band dispersion (top 
insets) . 



four Fermi points, for which we have a larger number of 
fc-points around Ep in a finite system. When [/ ^ A, 
dominance of the SDW is recovered as in the case of the 
ID chain and trestle lattice with smaller \t'\. 

Next we plot the exponents of the correlation functions 
as functions of t' in Fig. [TH The vertical dashed line in 
the figure represents the boundary at which the number 
of Fermi points changes from two to four. In other words, 
a geometrical frustration (i.e., interference between the 
nearest-neighbor and second-neighbor transfers) becomes 



the strongest around this boundary. We can notice that 
the pairing correlations (sSC and dSC) tend to be dom- 
inant around the boundary, while the CDW correlation 
begins to decay faster there. 



VII. SUMMARY 

To summarize we have shown for the ID Hubbard- 
Holstein model the following: (i) For the half-filled case 
we have obtained the phase diagram for the whole pa- 
rameter space spanned by the Hubbard U, phonon fre- 
quency Wo, and the electron-phonon coupling A. A re- 
gion is shown to exist between the SDW and CO phases, 
where the superconducting correlation is only subdomi- 
nant against CDW. (ii) When we either (a) dope the elec- 
tronic band, or (b) change the electronic band structure 
by considering a trestle lattice, the on-site pair correla- 
tion, and in the case of (b) the nearest-site singlet pair 
correlation, indeed become dominant over CDW. This is 
to be contrasted with the Tomonaga-Luttinger picture 
(g-ology), which would dictate that the pair correlation 
should not dominate when there are two Fermi points 
(the region left of the dashed line in Fig. [T4|) . 
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